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SYSTEMS AND METHODS FOR 
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WAVES USING CORRECTIVE 

SIGNAL 

Background of the Invention 

[0001 ] The systems and methods of the invention relate to the suppression of acoustic 
pressure waves, and in particular, to the suppression of acoustic pressure waves in 
gas turbine combustion chambers. 

[0002] It should be appreciated that adverse acoustic pressure waves may be generated 
in a variety of operational systems. For example, such adverse acoustic pressure wave 
may be generated in gas turbines. This problem may manifest itself when trying to 
increase the efficiency of the flames in the combustion chambers of such a gas 
turbine, for example. That is, the problem may manifest itself when trying to reduce 
the undesirable emissions generated by a gas turbine. By reducing the emissions, and 
the rate at which governmental emissions allotments are consumed, it is possible to 
maximize the number of revenue hours of the gas turbine. 

[0003] 

That is, when the flames are "leaned out," the emissions go down. However, the 
flame burning in the gas turbine may become unstable. Such an unstable flame 
creates a pressure wave, which may be of audible frequencies, and hence termed an 
"acoustic pressure wave." The acoustic pressure waves may stress various portions of 
the gas turbine, causing fatigue and shortening the life of the turbine. Specifically, 
torsional vibrations on the gas turbine's shaft may be created resulting in flexing and 
stressing the turbine blades. Additionally, the acoustic pressure wave may damage 
internal baffling in the combustion chamber. The acoustic pressure wave may also 
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adversely affect the efficiency of the machine. 

[0004] There are Icnown techniques relating to the active suppression of combustion 

chamber acoustics. However, the known techniques fail to teach an effective process 
to establish a corrective modulation signal at the correct magnitude, frequency and 
phase. Some known techniques create a modulation signal using adaptive filtering. 
The difficulty with such a technique lies in the need to, and time required for, filter 
coefficients to adapt when the acoustic signature is changing spectral content rapidly. 
Also, in the absence of any automatic gain control, the known techniques may actually 
exacerbate the undesired acoustic while re-adapting to the new spectral content. 
Accordingly, the known techniques suffer from the above drawbacks, as well as 
others. 

Brief Summary of the Invention 

[0005] The systems and methods of the invention solve the above problems, as well as 
other problems, present in known techniques. In accordance with one aspect, the 
invention provides a method for providing a corrective modulation signal to suppress 
an acoustic pressure wave in an operational system, the method comprising the steps 
of sampling the acoustic pressure wave generated in the operational system; sampling 
a previously generated corrective modulation signal, the previously generated 
corrective modulation signal having parameters; performing fast Fourier transform 
processing on the sampled acoustic pressure wave; performing a pair of single 
frequency discrete Fourier transform processing on the sampled acoustic pressure 
wave; determining the frequency, phase and magnitude of a dominate pressure wave 
in the acoustic pressure wave based on the fast Fourier transform processing and the 
discrete Fourier transform processing; and generating a sinusoidal corrective 
modulation signal to suppress the acoustic pressure wave based on the frequency, 
phase and magnitude of the dominate pressure wave and the parameters of the 
previously generated corrective modulation signal, the corrective modulation signal 
being at substantially the same frequency as, and generally 1 80 degrees out of phase 
with, the acoustic pressure wave. 

[0006] 1^ accordance with a further aspect, the invention provides a corrective 

modulation system for providing a corrective modulation signal to suppress an 
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acoustic pressure wave in an operational system, the system comprising a pressure 
sampling device that samples the acoustic pressure wave generated in the operational 
system to provide a sample of the acoustic pressure wave; a phase output portion, the 
phase output portion providing a sample of a phase of a previously generated 
corrective modulation; a signal processing portion that processes the sample of the 
acoustic pressure wave and the sample of the phase of the previously generated 
corrective modulation, the signal processing portion including a fast Fourier transform 
processing portion that performs a fast Fourier transform process on the sample of 
the acoustic pressure wave, the signal processing portion generating frequency with 
maximum power information and maximum power information based on the fast 
Fourier transform process;at least two discrete Fourier transform processing portions 
that perform single frequency discrete Fourier transform processing, the signal 
processing portion generating pressure phase information based on the single 
frequency discrete Fourier transform processing, and a modulation phase processing 
portion, the modulation phase processing portion generating modulation phase 
information based on the sample of the phase of the previously generated corrective 
modulation. The system further includes a corrective modulation generator that 
generates a sinusoidal corrective modulation signal to suppress the acoustic pressure 
wave based on the frequency with maximum power information and maximum power 
information, the pressure phase information, and the modulation phase information, 
wherein the corrective modulation signal being at substantially the same frequency as, 
and generally 1 80 degrees out of phase with, the acoustic pressure wave. 

[0007] 

In accordance with a further aspect, the invention provides a system for providing 
a corrective modulation signal to suppress an acoustic pressure wave in an 
operational system, the system comprising means for sampling the acoustic pressure 
wave generated in the operational system; means for sampling a previously generated 
corrective modulation signal, the previously generated corrective modulation signal 
having parameters; means for performing fast Fourier transform processing on the 
sampled acoustic pressure wave; means for performing a pair of single frequency 
discrete Fourier transform processing on the sampled acoustic pressure wave; means 
for determining the frequency, phase and magnitude of a dominate pressure wave in 
the acoustic pressure wave based on the fast Fourier transform processing and the 
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discrete Fourier transform processing; and means for generating a sinusoidal 
corrective modulation signal to suppress the acoustic pressure wave based on the 
frequency, phase and magnitude of the dominate pressure wave and the parameters 
of the previously generated corrective modulation signal, the corrective modulation 
signal being at substantially the same frequency as, and generally 1 80 degrees out of 
phase with, the acoustic pressure wave. 

[0008] In accordance with a further aspect, the invention provides a method for providing 
a corrective modulation signal to suppress an acoustic pressure wave in a gas turbine 
system, the method comprising the steps of sampling the acoustic pressure wave 
generated in the gas turbine system; sampling a previously generated corrective 
modulation signal, the previously generated corrective modulation signal having 
parameters; performing fast Fourier transform processing on the sampled acoustic 
pressure wave; performing a pair of single frequency discrete Fourier transform 
processing on the sampled acoustic pressure wave; determining the frequency, phase 
and magnitude of a dominate pressure wave in the acoustic pressure wave based on 
the fast Fourier transform processing and the discrete Fourier transform processing; 
generating a sinusoidal corrective modulation signal to suppress the acoustic pressure 
wave based on the frequency, phase and magnitude of the dominate pressure wave 
and the parameters of the previously generated corrective modulation signal, the 
corrective modulation signal being at substantially the same frequency as, and 
generally 1 80 degrees out of phase with, the acoustic pressure wave; generating a 
frequency error; generating a phase error; and providing a gain control based on the 
frequency error, the phase error and the magnitude of the dominate pressure wave, 
the gain control generating a gain signal to adjust the corrective modulation signal. 

[0009] 1^ accordance with a further aspect, the invention provides a corrective 

modulation system for providing a corrective modulation signal to suppress an 
acoustic pressure wave in an operational system, the system comprising a pressure 
sampling device that samples the acoustic pressure wave generated in the operational 
system to provide a sample of the acoustic pressure wave; a phase output portion, the 
phase output portion providing a sample of a previously generated corrective 
modulation; a signal processing portion that processes the sample of the acoustic 
pressure wave and the sample of the previously generated corrective modulation, the 
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signal processing portion including a fast Fourier transform processing portion that 

performs a fast Fourier transform process on the sample of the acoustic pressure 

wave, the signal processing portion generating frequency with maximum power 

information and maximum power information based on the fast Fourier transform 

process; at least one discrete Fourier transform processing portion that performs 

single frequency discrete Fourier transform processing, the single frequency discrete 

Fourier transform processing including performing a first single frequency discrete 

Fourier transform on a first part of the sample, which is processed by the signal 

processing portion to generate pressure phase information, and performing a 

second single frequency discrete Fourier transform on a second part of the sample, 

which is processed by the signal processing portion to generate pressure phase 

information, and a modulation phase processing portion, the modulation phase 

processing portion generating modulation phase information and modulation phase 

K 

information based on the sample of the previously generated corrective 

K— 1 

modulation; a corrective modulation generator that generates a sinusoidal corrective 
modulation signal to suppress the acoustic pressure wave based on the frequency 
with maximum power information and maximum power information; the pressure 
phase information and the pressure phase information; and modulation phase 
information and modulation phase information; and the corrective modulation 

K— 1 

signal being at substantially the same frequency as, and generally 1 80 degrees out of 
phase with, the acoustic pressure wave. 

[0010] 

In accordance with a further aspect, the invention provides a system for providing 
a corrective modulation signal to suppress an acoustic pressure wave in a gas turbine, 
the system comprising means for sampling the acoustic pressure wave generated in 
the gas turbine; means for sampling a previously generated corrective modulation 
signal, the previously generated corrective modulation signal having parameters; 
means for performing fast Fourier transform processing on the sampled acoustic 
pressure wave; means for performing a pair of single frequency discrete Fourier 
transform processing on the sampled acoustic pressure wave; means for determining 
the frequency, phase and magnitude of a dominate pressure wave in the acoustic 
pressure wave based on the fast Fourier transform processing and the discrete Fourier 
transform processing; means for generating a sinusoidal corrective modulation signal 
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to suppress the acoustic pressure wave based on the frequency, phase and magnitude 
of the dominate pressure wave and the parameters of the previously generated 
corrective modulation signal, the corrective modulation signal being at substantially 
the same frequency as, and generally 1 80 degrees out of phase with, the acoustic 
pressure wave; means for generating a frequency error; means for generating a phase 
error; and means for providing a gain control based on the frequency error, the phase 
error and the magnitude of the dominate pressure wave, the gain control generating a 
gain signal to adjust the corrective modulation signal. 

Brief Description of the Drawings 

[001 1] The present invention can be more fully understood by reading the following 

detailed description of the exemplary embodiments together with the accompanying 
drawings, in which like reference indicators are used to designate like elements, and 
in which: 

[001 2] Fig. 1 is a schematic diagram showing a system and process of generating a 

corrective modulation signal in accordance with one embodiment of the methods and 
systems of the invention; 

[001 3] Fig. 2 is a schematic diagram showing the details of the field programmable gate 
array, that generates the modulation waves, of Fig. 1 in further detail in accordance 
with one embodiment of the methods and systems of the invention; 

[001 4] Fig. 3 is a schematic diagram showing further details of the system of Fig. 1 , 

including the detailed implementation of the calculations shown in signal processing 
portion 20 of Fig. 1 in accordance with one embodiment of the methods and systems 
of the invention; 

[001 5] Fig. 4 is a diagram showing a comparison of the frequency range of an N element 
FFT bin versus the frequency range of a corresponding coarse single frequency DFT 
bin in accordance with one embodiment of the invention; 

[001 6] Fig. 5 is a diagram showing aspects of the phase frequency relationship of either a 
FFT or DFT bin in accordance with one embodiment of the methods and systems of 
the invention; 
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[001 7] Fig. 6 is a scliematic diagram showing a proportional integral control processing 
portion of Fig, 1 in further detail in accordance with one embodiment of the methods 
and systems of the invention; and 

[001 8] Fig. 7 is a diagram showing an acoustic pressure wave with harmonic distortion 
and a corrective modulation to be applied to suppress the acoustic pressure wave in 
accordance with one embodiment of the methods and systems of the invention. 

Detailed Description of tlie Invention 

[001 9] The systems and methods of the invention offer a technique for providing a 

corrective modulation signal to suppress an acoustic pressure wave in an operational 
system, such as a gas turbine, for example. In accordance with one embodiment of 
the invention, the method includes the steps of sampling the acoustic pressure wave 
generated in the operational system, performing a fast Fourier transform (FFT) on the 
sampled acoustic pressure wave, and performing two single frequency discrete Fourier 
transforms (DFTs) on the sampled acoustic pressure wave. 

[0020] The method further includes determining the frequency, magnitude and phase of 
the dominate spectral component of the acoustic pressure wave based on the FFT and 
the DFTs processing. Further, the method includes generating a sinusoidal corrective 
modulation signal to suppress the acoustic pressure wave at the same frequency and 
resulting magnitude as that of the dominate pressure wave. The phase of the 
corrective modulation signal is sampled and controlled in such a manner as to 
establish a 1 80 degree phase relationship, i.e., appropriately taking into account, 
propagation delay corrections, with the dominate spectral component of the acoustic 
pressure wave. 

[0021] Hereinafter, various aspects of the invention will be described in further detail. 

The systems and methods of the invention provide a corrective modulation signal for 
use in the suppression of acoustic pressure waves in an operational system, and in 
particular in a combustion chamber of a gas turbine. However, it should be 
appreciated that the invention is not limited to such application. That is, the method 
of the invention may be utilized in a variety of operating environments in which 
control of acoustic pressure waves is desired. 
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[0022] In accordance with embodiments of the methods and systems of the invention, a 
modulation is generated at the correct frequency and phase via a novel technique. 
This technique combines the spectral analysis of a Fast Fourier Transform (FFT) with 
the inherent phase information of a voltage-controiied oscillator implemented in a 
field programmable gate array. It should be appreciated that the method of the 
invention eliminates required hardware and reduces associated costs. Computational 
loading of the processing unit is also reduced, releasing this resource for other uses. 
Additionally, the spectral analysis of the pressure waves is made available to the 
turbine control system for protective actions, time tagging, trending, or further 
analysis, for example. 

[0023] Gas turbine combustion systems can experience dynamic pressure oscillations in 
the audible frequency range, i.e., acoustics. These oscillations, if of sufficient 
magnitude and persisting long enough, can damage the combustion system, reduce 
the life expectancy of the system and/or affect the operation of the turbine. It should 
be appreciated that active control of these acoustics poses certain signal processing 
problems. First, the dominant pressure wave must be identified both in terms of 
frequency and magnitude. This must be discerned from potentially, a very noisy 
spectral background. 

[0024] Secondly, a corrective modulation signal must be created for application to a 

secondary fuel valve or an air bleed valve, for example, or some other device used to 
actively affect turbine parameters in order to reduce or eliminate the acoustics. This 
corrective modulation signal is locked at a 1 80 degree phase relationship to the 
acoustic pressure wave that is be squelched, in other words, the corrective modulation 
signal is located so as to be generally 1 80 degrees out of phase with the acoustic 
pressure wave, i.e., so that the relationship serves to cancel out the undesired 
acoustic pressure wave. Additionally, it is desirable to adjust the magnitude of the 
modulation in direct proportion to the amplitude of the acoustic pressure wave and 
inversely proportional to the phase relationship of the modulation to the acoustic 
pressure wave. The adjustment of the magnitude of the modulation in direct 
proportion to the amplitude of the acoustic pressure wave allows the suppression 
efforts to be of a magnitude to reduce the acoustic, but no larger, thereby avoiding 
the introduction of other undesirable effects. The adjustment of the magnitude of the 
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modulation in inverse proportion to the phase relationship of the modulation to the 
acoustic pressure wave assures that until the desired 1 80 degree relationship is 
approximately established, the magnitude of the modulation will be small or zero. 
This avoids the undesirable case where, while locking, the modulation and acoustic 
are in phase and therefore the modulation is exacerbating the problem, 

[0025] The systems and methods of the invention provide for two primary objectives, 

which are met efficiently and in a real time manner. With reference to Fig. 1 , the first 
objective is to provide spectral analysis of an input signal coming from a differential 
pressure transducer 1 0 that monitors the pressure in a gas turbine combustion 
chamber, for example, in accordance with one embodiment of the methods and 
systems of the invention. However, it should be appreciated that other devices may be 
used to perform such input signal. In accordance with further aspects of the methods 
and systems of the invention, the analysis of the invention determines the spectral 
component of the acoustic pressure wave with the largest power, as well as the 
frequency associated with the acoustic pressure wave. The entire results, or a portion 
of the results, of the spectral analysis can also be made available to the rest of a 
turbine control system for protection or trending, for example. 

[0026] Secondly, the systems and methods of the invention provide for the generation of 
a sinusoidal modulation signal that is at the same frequency as the spectral 
component of the acoustic pressure wave with the largest power, and which is locked 
at a 1 80-degree phase relationship with that acoustic pressure wave. The systems and 
methods of the invention provide for meeting these objectives in a manner minimizing 
hardware and associated costs, and minimizing computational time. 

[0027] 

Fig. 1 is a block diagram showing aspects of the process in accordance with one 
embodiment of the methods and systems of the invention. A pressure transducer 10, 
which may be a differential pressure transducer, for example, is placed within the 
combustion chamber of a gas turbine, in accordance with one embodiment of the 
methods and systems of the invention. It should be appreciated that a transducer may 
be used that is not differential. Should the pressure transducer not be differential, the 
process of the invention can still be utilized by removing the D,C. part, i.e. the time 
invariant part of the signal by ignoring this component in the spectral analysis in a 
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manner described below. The output of the pressure transducer is passed through 
signal conditianing circuitry and then to an A/D converter 1 2, as shown in Fig. 1 . 

[0028] It should be appreciated that the frequency response range of the transducer may 
exceed the frequency range of interest for the acoustic pressure wave desired to be 
controlled. If the frequency response range of the transducer exceeds the frequency 
range of interest for the acoustic pressure wave, an anti-aliasing filter may be utilized 
in conjunction with the signal conditioning circuitry, i.e., before the A/D converter 1 2, 
Similarly, should significant broadband noise be present on the pressure transducer 
signal, an anti-aliasing filter may be required. 

[0029] The signal from the A/D converter 1 2 is simultaneously sampled along with a 

signal from a phase register 56, which travels along path 57, as shown in Fig. ] . The 
count of the phase register 56 gives the instantaneous phase of the current corrective 
modulation, i.e., the corrective modulation currently being applied and generated, in 
accordance with one embodiment of the methods and systems of the invention. This 
register is a part of a sinusoidal voltage-controlled oscillator (VCO), which generates 
the corrective modulation. The VCO is implemented in a field programmable gate 
array (FPGA) 50, the details of which are further shown in Fig. 2 and described below. 
That is, Fig. 2 is a schematic diagram showing the details of the field programmable 
gate array 50 that generates the modulation in accordance with one embodiment of 
the method of the invention. Accordingly, the field programmable gate array (FPGA) 
50 may be characterized as a corrective modulation generator 50. 

[0030] The simultaneous sampling of the A/D converter 1 2 and the instantaneous phase 
register 56 is done via a Direct Memory Addressing unit (DMA) 1 4 of a 
microprocessor. In accordance with one embodiment of the methods and systems of 
the invention, a total of 2048 pairs of simultaneous samples of the A/D converter 1 2 
and phase register 56 are taken, which provide input 1 6 and input 1 7 into the signal 
processing portion 20. However, it should be appreciated that the systems and 
methods are not limited to such sampling, i.e., variations of the 2048 sample 
arrangement may be utilized. The DMA 14 allows the sampling to proceed without any 
participation by a main or a central processing unit, for example. Therefore, sampling 
can occur in parallel with the processor computing values on a last set of samples. 
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[0031] In accordance with one embodiment of the methods and systems of the invention, 
firmware, running on a micro processor, processes the 2048 pairs of samples. The 
first steps of such processing are shown in the signal processing portion 20, shown in 
Fig. 1 . Also, the process is shown in further detail in Fig. 3 and described below. The 
signal processing portion 20 includes the processing portion 22, the processing 
portion 23, the processing portion 24, the processing portion 25 and the processing 
portion 26, as shown in Fig. 1 . 

[0032] In accordance with one embodiment of the methods and systems of the invention, 
a mathematical windowing algorithm is used on the 2048 samples output from the 
A/D converter 1 2. This windowing is necessary to prevent signal discontinuities at the 
beginning and the end of the sampling from being analyzed as high frequency 
components. Such components, along with possible aliasing, might appear as false 
spectral components anywhere in the analyzed spectrum. The windowing serves to 
shape the samples, forcing the samples to zero at the first and last sample. As a 
result, discontinuities at the ends of the sampling may be removed. In accordance 
with embodiments of the methods and systems of the invention, a windowing process 
may be performed using one of a variety of known methods including the Rectangular 
method, the Hamming method, the Manning method, the Triangular method, the 
Blackman method, the Blackman - Harris method or the Flat Top method, for 
example. 

[0033] The output from the windowing process, in accordance with one embodiment of 
the methods and systems of the invention, then has a Fast Fourier Transform (FFT) 
performed on such output, as illustrated in processing portion 22 in Fig. 1. Next, for 
each complex element in the FFT, the power is calculated in processing portion 23, as 
shown in Fig. 1 , according to the equation: 

POWER = [(REALPARTOF FFT)' + (IMAGINARY PART OF FFT)' J/ FFT LENGTH' 

(Equation 1) 

[0034] 1^ accordance with one embodiment of the methods and systems of the invention, 
the maximum power of all the FFT elements is determined and referred to in the 
processing portion 23 as MAX POWER. The frequency associated with this power, 
which is referred to as FREQ WITH MAX POWER in processing portion 23 in Fig. 1 , is 
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also known, i.e., since each element E, in the FFT, has a frequency associated with it 
according to the equation: 

FREQUENCY ^ = ~ * SAMPLING FREQUENCY 

^ FFT LENGTH 

(Equation 2) 

[0035] where: 

[0036] E = FFT element number, also referred to as a bin number, and ranges from 0 to 
the (FFT length - 1), which in accordance with one embodiment of the invention is 0 
to 2047; and 

[0037] FFT LENGTH - number of samples on which the FFT is performed, which in 
accordance with one embodiment of the invention is 2048 , 

[0038] It should be appreciated that since the input is real and not complex, only 

elements E = 0 to E = (FFT length / 2) are independent and therefore the FFT bin 
frequencies should be computed over such range. Therefore, in accordance with one 
embodiment of the methods and systems of the invention, the FFT bin frequencies will 
range from 0 to a Sampling Frequency/2, and associated bin numbers from 0 to 1024. 

[0039] As noted above, it should be appreciated that rather than the differential pressure 
transducer 10, a pressure transducer may be utilized that is not differential. If the 
pressure transducer is not differential, the first FFT bin can be eliminated from the 
maximum power determination. This bin contains the steady or D.C. component and 
is generally not of concern in acoustic suppression, in accordance with embodiments 
of the methods and systems of the invention. 



[0040] 



It should be appreciated that the attenuation of the magnitude near a frequency 
bin's boundaries or edges is affected by the windowing selected. This "roll off," 
however, Is typically present to some degree and is pictured in Fig. 4. Fig. 4 is a 
diagram showing aspects of a single frequency bin in accordance with one 
embodiment of the invention. Further, Fig. 5 illustrates the phase shift that must be 
corrected across a frequency bin of any Fourier transform. Mathematically the phase 
shift occurring across a frequency bin may be calculated as follows: 
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PHASE SHIFTIN DEGREES = 180* ( K- M ) * ( 1 - ( 1 /N ) ) 



(Equation 3) 



[0041] 



Where: 



[0042] 



M = bin number, with bin numbers starting at 0; 



[0043] 



K = number of cycles in n samples of the frequency of interest; and 



[0044] 



N = number of samples. 



[0045] It should be appreciated that both of these phenomena cause difficulty in 

accurately calculating the phase of a frequency component that is at or near the edges 
of a frequency bin. To alleviate this situation, two single frequency DFTs are 
performed as shown in processing portion 24 of Fig. 1 . Each single frequency DFT is 
calculated at the frequency of the maximum power that has already been determined. 
Each single frequency DFT is also performed on exactly half the samples that the FFT 
was performed on, i.e., 1024, in accordance with one embodiment of the methods 
and systems of the invention. This will result in the width of the single frequency 
DFTs frequency bin being twice as wide in terms of frequency as the frequency bins 
of the FFT. The resulting resolution in the frequency spectrum is illustrated in Fig. 4. 
This shows a simple example of the mapping of the full length FFT into the more 
coarse frequency resolution of the single frequency DFT. The bin resolution illustrated 
in Fig. 4 is not necessarily the preferred frequency spectrum bin resolution, but is 
merely illustrative, in accordance with one embodiment of the methods and systems 
of the invention. As shown in Fig. 4 there is a direct mapping of the even frequency 
bins of the full length FFT into the corresponding single frequency DFT bin with the 
relationship: 

Single Frequency DFT Bin Number = Full Length FFT Bin Number / 2 

(Equation 4) 

[0046] 

However there exists an ambiguity with the mapping of the odd full length FFT 
frequency bins into the single frequency DFT frequency bins as shown in Fig. 4 by the 
thick horizontal bans on either side of the single frequency DFT bin boundaries. As 
can be seen by the dual arrows mapping the full length FFT bin to either side of the 
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single frequency DFT bin boundary a choice has to be made when mapping the odd 
numbered full length FFT bins: 

Single Frequency DFT Bin Number = truncate (Full Length FFT Bin Number 12) 



[0047] By default, Equation 5 is the initial attempt to map the odd full length FFT 

frequency bin into the single frequency DFT frequency bin, in accordance with one 
embodiment of the methods and systems of the invention. If a frequency lock is not 
achieved within a predetermined number of scans the alternate mapping, Equation 6, 
is implemented to acquire lock. 



[0048] Hereinafter, further aspects of the 2048 samples will be described in accordance 
with one embodiment of the methods and systems of the invention, and in particular, 
processing relating to the single frequency DFTs. The 2048 samples may be 
considered as two groups as shown in the processing portion 24 of Fig. 1 . The oldest 
1 024 samples, i.e. the first 1 024 samples taken, may be characterized as group K-1 . 
The newest 1024 samples, i.e. the last 1024 samples taken, may be characterized as 



group K. Each group of samples is windowed in a manner similar to that done on the 
entire 2048 samples as discussed above. Each windowed group of samples then has a 
single frequency Discrete Fourier Transform performed on it at the frequency with the 
maximum power that has been calculated, as discussed above. This processing is 
shown in the processing portion 24 in Fig. 1 , as well as Fig. 3, and is designated as 
"SINGLE FREQ DFT." 

[0049] Further, the phase angle of the of the frequency at which the maximum power was 
found is calculated, in processing portion 25, for group K as follows: 



(Equation 5) 



or 



Single Frequency DFT Bin Number = (Full Length FFT Bin Number / 2) + 1 



(Equation 6) 




(Equation 7) 



[0050] 



The meaning of this phase angle should be appreciated. That is, such phase angle 
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is the phase of the spectral component, who's frequency at which the maxinnum 
power was found, at the instant that the first sample of FFT was taken. 

[0051] In further description of the systems and methods of one embodiment of the 

invention, the phase angle of the frequency at which the maximum power was found, 

at the instant that the first sample of DFT was taken, may now be calculated. This 

K— 1 

calculation is performed, by the processing portion 25, using group K-l. The equation 
for this processing in accordance with one embodiment of the methods and systems 
of the invention is: 

( IMAGINARY PART OF SINGLE FREQ DFTk., 



PHASE ^ ARC TAN 



REAL PART OF SINGLE FREQ DFT^., 

(Equation 8) 

[0052] It should be appreciated that if there is significant propagation delay from when 
the modulation is applied to when an effective change in the acoustic results, a 
compensation can be made. This compensation is shown in Fig. 3, i.e., processing 
portion 322, and is some number of degrees being subtracted from the modulation 
phases. The compensation degrees can be calculated as follows: 

[0053] Propagation delay compensation in degrees = 

VCO FREQ IN HZ. * 360 DEG. / CYCLE * Propagation delay in seconds 

(Equation 9) 

[0054] Now, remembering that the A/D converter samples and those of the instantaneous 
phase register were taken simultaneously, it becomes easy to calculate the phases of 
the corrective modulation that correspond to the PHASE K and PHASE K-l , which were 
just calculated. That is, one merely selects sample number 1 and sample number 
1025 from the instantaneous phase register and, thereafter, scales the count values 
into degrees. These scaled values, with compensation for propagation delays if 
needed, become the modulation phase angles referred to in Fig. 1, i.e., in processing 
portion 26, and Fig. 3, as MODULATION PHASE and MODULATION PHASE . 



[0055] 



With further reference to Fig. 1 , two phase errors, PHASE ERROR 28 and PHASE 

K 

ERROR 27, as shown in Fig. 1 , may then be calculated by merely subtracting the 
K— 1 
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A/D phases from the corresponding modulation phases. Any necessary phase 

correction, due to the difference between the VCO"s frequency and the center 

frequency of the associated FFT bin, is calculated at processing portion 30 in Fig. 1 , 

PHASE ERROR is then corrected in calculation portion 31 for the 90 degree shift 
K 

between the VCO sine based phase and the cosine based phase, which are calculated 
by the single frequency DFT, 

[0056] Additionally, the PHASE ERROR is corrected in the calculation portion 31 for the 
phase that occurred due to the frequency component's placement within the 
frequency bin. Finally, the PHASE ERROR is shifted 1 80 degrees in calculation 
portion 32, fed into a proportional and integral control 33 with its gains of G (Phase 
path. Integral gain) and G (Phase path, Proportional gain). As a result, the method 
of the invention produces the CORRECTED PHASE ERROR 34, as shown in Fig. 1 , The 
schematic diagram of Fig. 6 details the proportional and integral control (33, 39) used 
both in the phase and the frequency paths of the invention. 

[0057] In accordance with embodiments of the methods and systems of the invention, the 
difference of the two phase errors i.e., PHASE ERROR ^ PHASE ERROR is compared 
in calculation portion 36 and is used to calculate a slip frequency, or what may be 
characterized as a FREQUENCY ERROR 35 in Fig. 1 . It should be appreciated that since 
the change in phase error occurred over the time it took to accumulate the number of 
samples for a DFT, the difference is first multiplied by 1 over that time as illustrated in 
calculation portion 37 of Fig. 1 , i.e. 

[0058] ^/T^^ 

[0059] Next, in accordance with one embodiment of the methods and systems of the 
invention, in order to change from units of degrees/time to cycles/time a 
multiplication by 1 /360 is done in calculation portion 38 of Fig. 1 . The result is fed 
into a proportional and integral control 39, as well as the gains of C (Frequency 



[0060] 



path. Integral gain) and G (Frequency path, Proportional gain). 



The output of the proportional and integral control 39 is added in calculation 
portion 40 to the count value residing in a counter 42, i.e., which counts the number 
of times the phase integrator hit either a positive or negative clamp. In accordance 
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with one embodiment of the methods and systems of the invention, a count of (+1) is 
added to the counter 42, by a processing portion 41 as shown in Fig. 1 , every time the 
phase integrator hit the positive clamp and the integrator was reset, and a count of (- 
1 ) is added for each time the negative clamp was hit and the Integrator was reset. The 
sum of the counter 42 and the output of the frequency proportional and integral 
control 39 forms the corrected frequency error 35, as shown in Fig. 1 . 

[0061] It should be appreciated that it now becomes possible to determine the 

instantaneous frequency at which the VCO is to run. That Is, the sum of the corrected 
frequency error 35 and the corrected phase error 34 is subtracted from the FREQ 
WITH MAX POWER output and then scaled into the appropriate count value to make 
the VCO run at that frequency, i.e., the V.C.O. frequency 55 as shown in Fig. 1 . In 
accordance with one embodiment of the methods and systems of the invention, the 
equation for the VCO frequency with the 1 MHZ clock 206 and register 204, as shown 
in Fig. 2, and accumulator widths, as shown in Fig. 2, is as follows: 
/ 



16384 



(l * 10"^ SEC)2^^^ ^ EXECUTION CLOCK DIV REG) 



1 

* 1024 



CNT IN FREQ SELECTION REG 

(Equation 10) 

[0062] Solving this equation for the counts to be placed in the frequency selection 
register results in the relationship: 

CNT IN FREQ SELECTION REG = 

1 6 77721 6 * 2 execution cux:k div reg ) * p^^^ 

(Equation 11) 

[0063] It should further be appreciated that the capability to divide down the 1 MHZ clock 
206 may also be provided using a suitable input 202, as shown in Fig. 2. Alternatively, 
the execution clock divide register 204 may simply be fed a constant one. 



[0064] 



Also, having caused the VCO to run at the correct frequency for locking the 
modulation, the amplitude of the modulation needs to be addressed. It is desirable to 
increase the amplitude of the modulation when the acoustic pressure wave increases 
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in magnitude. It is also desirable to have little or no amplitude if tlie modulation is not 
dose to a 1 80 phase shift with respect to the acoustic pressure wave. This should be 
appreciated by one of ordinary skill in the art, since otherwise the modulation will only 
worsen the acoustic. In accordance with one embodiment of the methods and systems 
of the invention, the solution to both these desired actions is to add an Automatic 
Gain Control (ACQ 46, as shown in Fig. 1 . 

[0065] By providing the AGC 46 with both the maximum power of the acoustic pressure 
wave 70, the frequency error 72, and the phase error 74, the AGC 46 can adjust the 
amplitude of the modulation as desired. Specifically, the AGC 46 outputs a gain signal 
via the path 53 to the FPGA 50. The gain signal 53 is proportional to the maximum 
power of the acoustic pressure wave and inversely proportional to the absolute value 
of the phase error and the absolute value of the frequency error. Priority is given to 
the phase error so that no modulation is provided until the phase relation nears 1 80 
degrees regardless of the magnitude of the acoustic pressure wave. 

[0066] Hereinafter, further features of the systems and methods of the invention will be 
described with further reference to Fig. 2 and implementation of the VCO in the FPGA 
50. It should be appreciated that a core component of the VCO is a 24 bit wide 
accumulator register 216, as shown in Fig. 2. An accumulation execution occurs at a 

(CNT IN EXECUTION CLOCK DIV REG) . 
rate equal to the 1 MHZ clock 206 divided by 2 m 

accordance with one embodiment of the methods and systems of the invention. At 
this time the contents of the accumulator 21 6 are added to the contents of a 
frequency selection register 21 0 and the sum is then output to accumulator 216. Bits 
14-23 of the accumulator 21 6 are then mapped to bits 0-9 of the instantaneous 
phase register 56. Accumulator bits 14 22 are used as an index into a sine magnitude 
table 21 2, as shown in Fig. 2. The table 212 contains the magnitudes of a sine wave 
for the range of 0 to 1 79.64844375 degrees, in accordance with one embodiment of 
the methods and systems of the invention. 

[0067] 

In accordance with one embodiment of the invention, the contents of the table 
21 2 are counts from 0 to 255 corresponding to 0 to 1 .0. The value of the indexed 
table entry, i.e., bits 0 - 8 is mapped as bits 2 - 1 0 of a shift register 21 8. Bits 0 and 1 
of the shift register 21 8 are set to 0. The half sine wave contained in the table 21 2 is 
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expanded to a full sine wave by the use of bit 23 of the accumulator 216, which maps 
to bit 1 1 of the shift register 21 8. Using this technique reduces real estate, i.e., 
memory , utilization for the sine table 21 2 in the FPGA 50. Finally, the contents of the 
shift register 21 8 are shifted right by the count value in the magnitude selection 
register 54 and passed to the D/A 60 via path 62. 

[0068] Hereinafter, further aspects of the systems and methods of the invention will be 
described with reference to Fig. 3. Fig. 3 is a diagram showing further details of the 
processing being performed in the signal processing portion 20 of Fig. 1 . As shown in 
Fig. 3, input 302 includes the two samples that are simultaneously taken, by the DMA 
14, from the A/D 1 2 and from the F.P.G.A.'s Instantaneous phase register 56. These 
two samples from the A/D 1 2 and from the phase register 56 are input into the 
memory portion 306 and the memory portion 316, respectively. Each memory portion 
306 and 31 6 is provided with 2048 respective samples. 

[0069] The samples in memory portion 306 are then passed through respective windows. 
Specifically, the window 308 is used in conjunction with the calculation of the (K-1) 
single frequency DFT in the processing portion 312. The window 31 0 is used in 
conjunction with the calculation of the (K) single frequency DFT in the processing 
portion 314. Also, the window 328 is used in conjunction with the calculation of the 
FFT in the processing portion 330. 

[0070] The windows 308, 31 0 and 328 utilize a mathematical windowing technique, such 
as a Blackman or Flat Top technique, as is described above. It should be appreciated 
that the windows 308, 31 0 and 328 may use the same windowing technique or 
different windowing techniques. 

[0071 ] In accordance with this embodiment of the invention, further aspects of the FFT 
processing will be described. As noted above, a complete set of 2,048 samples, i.e., 
2048 samples of the acoustic pressure wave at 2048 addresses, is available in the 
memory 306. The desired output of the window 328 and the FFT processed in the 
processing portion 330 is a power versus frequency matrix, and a phase versus 
frequency matrix or relationships. Specifically, the pressure signal is broken into 
spectral or frequency components using the window 328, the processing portion 330 
and the processing portion 332, aspects of which are also described above. 



Page 19 of 41 



[0072] Each of these spectral components is then looked at by the processing portion 
334 to determine the magnitude at each of the frequency components. Then, the 
processing portion 334 determines the frequency component with the largest power 
and therefore, the largest magnitude. The relationship of magnitude to power is given 
by the following equation: 

[0073] MAGNITUDE = SQUARE ROOT ( POWER ) 

[0074] As a result, the processing portion 334 generates output 338 providing the 

maximum power and the output 340 providing the particular frequency that has that 
maximum power. 

[0075] As a result, the power, and magnitude of the acoustic is known and the 

approximate frequency of the acoustic is known. However, it should be appreciated 
that the FFT processing produces frequencies with only a certain resolution. 
Accordingly, the frequency of the acoustic is only approximately known and the power 
and magnitude are known. Thus, it should be noted that the frequency is not known 
exactly, nor has phase information been determined. These further determinations are 
provided by the other processing of Fig. 3. 

[0076] Hereinafter, operations of the processing portion 31 2 and the processing portion 
314 will be further described. In accordance with this embodiment of the invention, 
the 2048 samples are split into two sets of 1024. A windowing process is performed 
on each set using the windows 308 and 31 0. Further, a single frequency DFT is 
performed on each set in the processing portion 31 2 and the processing portion 314, 
respectively. The single frequency DFT"s are calculated for the frequency at which the 
maximum power occurred as determined by the FFT , i.e. the frequency with 
maximum power 340, as shown in Fig. 3. This processing yields a phase. To explain 
further, the window 308 and the processing portion 312 yield an instantaneous phase 
of the single frequency at the time of the first sample, i.e., sample one, out of the 
2048 list. Further, the window 310 and the processing portion 314 yield an 
instantaneous phase of the single frequency at the time of sample 1 025. Accordingly, 
the acoustic"s phase is measured and thus known at two points in time as represented 
by sample 1 and sample 1 025. 
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[0077] The output from the processing portion 31 2 is input into the processing portion 
324. Also, the frequency with the maximum power information 340, described above, 
is input into the processing portion 324. Based on this input, the processing portion 
324 determines the phase of the frequency with the maximum power based on the 
first set of samples, and outputs this information as output 335, The output 335, 
which is based on the first sample, thus may be characterized as the acoustic pressure 
phase (K-1). 

[0078] Also, the output from the processing portion 314 is input into the processing 

portion 326. Also, the frequency with the maximum power information 340, described 
above, is input into the processing portion 326. Based on this input, the processing 
portion 326 determines the phase of the frequency with the maximum power based 
on the second set of samples, and outputs this information as output 336. The output 
336, which Is based on the second sample, thus may be characterized as the acoustic 
pressure phase (K). 

[0079] Turning now to the samples in the memory 316, which are input from the phase 
register 56, these samples provide the instantaneous phase angles of the 
modulations, i.e., including the instantaneous phase angles, for sample 1 and sample 
1025. These samples are taken simultaneously with the pressure waves. Thus, the 
method of the invention has allowed a phase angle comparison to be made. That is, at 
two points in time it is known what the acoustic phase was and what the modulation 
phase was. 

[0080] 

Accordingly, the process of the invention generates an output 318 based on 
sample 1025, which provides modulation phase K information. Also, the output 320 
based on sample 1 is generated that provides modulation phase (K-1) information. 
These outputs are adjusted using a scaling factor 317, i.e., "scaling in Deg/CNT." This 
is done since all the samples in the memory 316, of Fig. 3, are in units of counts. 
These counts have a specific conversion to units of degrees. In one exemplary 
implementation of the invention, 1 FPGA count corresponds to 0.351 5625 degrees. 
Also, the modulation phase outputs based on sample 1 and sample 1025 may be 
adjusted by a propagation delay compensation factor. This propagation delay 
compensation factor may be calculated by the processing portion 322, as is described 
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above. 

[0081] Fig. 5 is a diagram showing in further detail the proportional and integral control 
portion 33 of Fig. 1 , in accordance with one embodiment of the invention. As shown 
in Fig. 6, an input signal 602 is input into the proportional and integral control 
portion 33 as described above. This input signal is then split into two signals, i.e., 
signal 604 and signal 606. A gain G ^ 608 is applied to the signal 604, thus resulting 
in the adjusted signal 61 6. Further, a gain G ^ 61 0 is applied to the signal 606. These 
respective gains G p and G ^ may be determined in any suitable manner. As used "G ^ " 
is the symbolic name for "Gain of the Proporttonal" and "G ^ " is the symbolic name for 
"Gain of the Integral." 

[0082] As shown in Fig. 6, once the gain 61 0 is applied to the signal 606, then a suitable 
transfer function is applied in the processing portion 612. A transfer function defines 
the relationship between the inputs to a system and its outputs. The transfer function 
is typically written in the frequency, or V domain, rather than the time domain. 
Illustratively, a Laplace transform, for example, may be used to map the time domain 
representation into the frequency domain representation. The specific transfer 
function 1 /S represents an integration. Once the transfer function is applied, a clamp 
614 is then applied to the signal, thus resulting in the adjusted signal 61 8. The 
adjusted signal 61 8 is then added to the adjusted signal 61 6 in order to generate an 
output signal 620. The output signal 620 is then used in further processing in 
accordance with some embodiments of the invention, as is described above. 

[0083] Fig. 6, as described above, shows further details of the proportional and integral 
control portion 33. It should be appreciated that the proportional and integral control 
portion 39, as shown in Fig. 1 , may utilize a similar arrangement as that of the 
proportional and integral control portion 33. Accordingly, details of the proportional 
and integral control portion 39 have not been described in further detail. 

[0084] |p further explanation of the systems and methods of the invention. Fig. 7 is a 
scope tracing showing an acoustic pressure wave 402 with harmonic distortion, as 
well as a corrective modulation wave 404, which is generated using the process of the 
invention. Specifically, the acoustic pressure wave 402 is generated from the actual 
signal coming from the differential pressure transducer 1 0, as shown in Fig. 1 . It 
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should be appreciate that in accordance with the methods and systems of the 
invention, the corrective modulation wave 404 may be scaled in such a manner so as 
to essentially eliminate the acoustic pressure wave 402. That is, the two waves cancel 
each other out. It should be noted that, as seen in Fig. 7, the distortion at the troughs 
and crests of the respective waves is an artifact of a digital scope that was utilized, as 
opposed to an actual frequency component in the respective waves. 

[0085] It should further be appreciated that it may be desirable to not achieve 1 00% 

cancellation of the acoustic pressure wave, but rather to reduce the magnitude of the 
wave to a magnitude where it possesses insufficient power to cause damage or 
adversely affect machine performance. That is, in accordance with one embodiment of 
the methods and systems of the invention, a residue of the acoustic pressure wave is 
allowed to remain. This allows the corrective modulation 404 to remain "locked on" to 
the acoustic pressure wave 402. In contrast, if the acoustic pressure wave 402 was 
completely canceled out, the lock-on would be lost. Accordingly, it may be desired to 
control the magnitude of the corrective modulation. This control is performed using 
the automatic gain control 46, as described above. 

[0086] Also, the systems and methods in accordance with various embodiments of the 
invention have been described above using 2048 samples. However, as noted above, 
it should be appreciated that the practice of the invention is not limited to such a 
sample size. Rather, other suitable sample sizes may also be used. 

[0087] While the foregoing description includes many details and specificities, it is to be 
understood that these have been included for purposes of explanation only, and are 
not to be interpreted as limitations of the present invention. Many modifications to 
the embodiments described above can be made without departing from the spirit and 
scope of the invention, as is intended to be encompassed by the claims and their legal 
equivalents. 
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